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Abstract 

We introduce and investigate numerically a minimal class of dynamical tri- 
angulations of two-dimensional gravity on the sphere in which only vertices of 
order five, six or seven are permitted. We show firstly that this restriction 
of the local coordination number, or equivalently intrinsic scalar curvature, 
leaves intact the fractal structure characteristic of generic dynamically triangu- 
lated random surfaces. Furthermore the Ising model coupled to minimal two- 
dimensional gravity still possesses a continuous phase transition. The critical 
exponents of this transition correspond to the usual KPZ exponents associated 
with coupling a central charge c = ^ model to two-dimensional gravity. 



That two-dimensional {2D) quantum gravity can be regularized and studied 
using dynamical triangulations (DT) is well known Q and has lead to extensive 
study of the properties of dynamically triangulated random surfaces (DTRS). In- 
deed, many of the properties of critical spin systems on such lattices have been 
shown to follow from continuum treatments of central charge c < 1 theories coupled 
to 2D-quantum gravity. Furthermore, this formulation has the advantage that it 
renders tractable the calculation of many features of these theories which are inac- 
cessible to continuum methods. Foremost amongst these are questions related to 
the quantum geometry and fractal structure of the typical 2D manifolds appearing 
in the partition function [0. 

It has often been speculated that some of these models should also find a real- 
ization in condensed matter physics as models of membranes and/or interfaces with 
fluid in-plane degrees of freedom. The problem in trying to map such problems onto 
dynamically triangulated models, however, has always been the occurrence of large 
vertex coordination numbers in an arbitrary random graph ~ real condensed matter 
systems typically place a cut-off on the number of interactions allowed by the mi- 
croscopic degrees of freedom. It is possible that such cut-offs are unimportant in the 
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critical region of these theories. As a crude analogy we note that the short-distance 
differences between different discrete polygonal decompositions of two-dimensional 
gravity are, without fine-tuning, irrelevant in the continuum limit 

In this letter we present numerical results for a new model (called the minimal 
dynamical triangulation or MDT model) which is based on random triangulations 
but with an additional constraint: the coordination number of any vertex in the 
lattice can only be five, six or seven. Such minimal models implement a severe form 
of the microscopic interaction cut-off required by condensed matter systems and 
correspond in gravitational language to the inclusion of only three possible states 
for the local curvature — 27r/7, 0, 27r/5. (The scalar curvature at a vertex with 
coordination number qi is given by i? = 27r^^.) We have studied two systems: 
pure gravity as modelled by pure (minimal) triangulations and a set of Ising spins 
coupled to a dynamical lattice generated by minimal triangulations. We show in 
both cases that the critical behavior of the system appears to be identical to the 
conventional DTRS models - yet another striking indication of universality. 

The dynamically triangulated version of 2D-gravity is given by the partition 
function 

Z(^) = ^e-^^Z(A), (1) 

A 

Z{A) = Y.Zm. (2) 
TeT 

Z(A) is the fixed area partition function, fj, the cosmological constant and T the 
class of triangulations summed over. Zm represents the partition function of some 
generic matter field living on the surface. A given choice of T corresponds to a par- 
ticular discretization of the surfaces. Two classes commonly used are combinatorial 
and degenerate triangulations, Tc and T^. In the former triangles are glued to- 
gether along edges so as to form closed manifolds, with the constraints that any two 
triangles can only have one edge in common and no triangle can be glued onto itself. 
For degenerate triangulations these constraints are relaxed - this amounts to allow- 
ing tadpoles and self-energy diagrams in the dual representation to triangulations. 
Clearly TcCTd. 

A priori different classes T do not have to lead to the same theory, although 
the two aforementioned do, at least for the case of Ising matter where the models 
are exactly soluble |Q . Clearly this would no longer be true if T only contained one 
triangulation - an extreme limit. The class of triangulation Tm we investigate in 
this letter is a subclass of Tc defined by allowing only vertices of order 5, 6 or 7. 
This is as far one can go in suppressing curvature fluctuations on the surface while 
still retaining its dynamical (fluid) nature. 

We start by investigating the minimal model in the absence of matter {Zm = 1 
in Eq. 2). To see if this model belongs to the same universality class as pure 
gravity with triangulation space 7^ or To we study its critical behavior numerically. 
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Figure 1: The measured distributions of baby universes for (a) the minimal class of 
triangulations, (h) combinatorial triangulations and (c) triangulations obtained by 
node decimation from (a). All measurements are for triangulations of 1000 nodes. 

More precisely we measure the string susceptibility exponent 7^, defined through 
the scaling Z{^) ~ (/Xc — ^)'^~^" ■ The method we employ is that of calculating the 
distribution of minimal baby universes n{B) on the ensemble of triangulations Q. 
A baby universe of area B \s a. region of the triangulation joined to the bulk solely 
by a minimal neck. For the class Tm (and Tc) a minimal neck is a loop consisting 
of three links. Using the expected asymptotic behavior of the fixed area partition 
function, Z{A) ~ exp(;Uc^)^'^°~^) a simple argument relates n{B) to 7^ |§, 

nA{B) ~ [B{A - B)y'~^ , A ^00. (3) 

We have simulated the fixed area partition function Z{A), using Monte Carlo 
methods, for triangulations consisting of 1000 and 4000 vertices Q. Henceforth 
we will use the number of vertices as our measure of the area A. In Fig. 1 we 
show the measured distribution n{B) for A = 1000. For comparison we include 
the corresponding distribution measured on the class of combinatorial manifolds 
7(7 • Fig. 1 also includes a third distribution (c) whose significance will be addressed 
shortly. The first thing we notice is that these distributions are very different for 
B < 100 - in particular there are no small baby universes for the class Tm as the 
curvature restriction puts a lower limit on the size of baby universes; our construction 
effectively smoothens the triangulation locally. 

For larger baby universes the distributions look much more alike. To measure 7^ 
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Table 1: Measured values of 7^ for minimal dynamical triangulations (without mat- 
ter) after applying varying levels of node decimation with a blocking factor of 6 = 2. 
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we fit the tail of n{B) to Eq. 3. This gives 7, = -0.64(5) {A = 1000) and -0.53(3) 
{A = 4000). This agrees with the exact value for pure gravity 7^ = —0.5. The 
corresponding value obtained using 7^ is 7^ = —0.501(4) {A = 1000) and —0.504(4) 
[A = 4000). It is easy to understand why the values for minimal triangulations 
are less accurate - finite-size effects are more pronounced since more data for small 
values of B has to be excluded before we reach the asymptotic regime. Yet the result 
indicates that the minimal triangulated model is in the same universality class as 
the full-blown DTRS model. 

Additional evidence for this claim is obtained by applying a recently proposed 
Monte Carlo renormalization group method for blocking dynamical triangulations 
{node decimation ||8|) to the model. Each triangulation is blocked by removing 
vertices at random and in the process the restrictions on the curvature are dropped 
so that the model flows into the wider class of combinatorial triangulations. In this 
way we can ascertain whether the model after blocking is the same as that obtained 
with the triangulation class Tc - in particular we would like to know if the RG 
transformation generates a flow towards a fixed point, accessible from combinatorial 
triangulations, and corresponding to pure gravity. 

Curve (c) of Fig. 1 is the measured distribution n{B) for surfaces of 1000 vertices 
obtained from surfaces of 4000 vertices by node decimation. It is apparent that this 
distribution is much closer to the distribution for combinatorial triangulations than 
the minimal one. This we take as an indication that the model flows towards the 
same non-trivial pure gravity flxed point as combinatorial manifolds. This evidence 
can be strengthened by measuring the exponent 7^ on the set of triangulations 
obtained with blocking. This we show in Table || for various degrees of blocking. 
Agreement with the pure gravity value clearly improves with blocking. This conflrms 
that the model, after blocking, is closer to the pure gravity fixed point. 

The RG fiows above can be visualized by looking at the expectation values of 
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Figure 2: The RG flow of curvature operators in the (Oi, 02)~plaiie. The tra- 
jectory for minimal triangulations is marked with solid circles, while open symbols 
correspond to combinatorial manifolds with varying degrees of coupling to a loga- 
rithmic curvature term. The dashed circle indicates a non-trivial fixed point, while 
the solid one is a zero-volume sink. 

a subset of operators such as 



where qi is the curvature of a node i, k is the blocking level and the first sum is 
over nodes while the second is over nearest neighbors. Different trajectories in the 
(Oi, 02)-plane can be accessed by adding an irrelevant operator to the bare action 
and varying its coupling constant. We have chosen aJ2i^*^SQi ■ This is shown in 
Fig. 2. By changing the coupling constant a we effectively change the trajectories 
identifying the location of the non-trivial fixed point. And indeed the trajectory for 
the minimal triangulations flows directly towards the pure gravity fixed point before 
it veers off and follows a universal trajectory into a zero-volume sink. 

Our findings thus far are consistent with the exact solution of two-dimensional 
gravity given in Q. These authors solve a matrix model of dually weighted 
planar graphs which allows the incorporation of the higher-curvature operator R^. 
They find that the R^ operator is non-perturbatively irrelevant in the continuum 
limit. The model always reduces at large length scales to a model of pure gravity 




and 
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- there is no transition to a "flat" phase. For technical reasons the solution of |^ 
is restricted to lattices with even coordination number vertices. There is thus no 
direct mapping to the model we treat. Nevertheless in both cases it is found that 
higher-curvature contributions to the effective action are irrelevant. 

In the second part of this letter we investigate the Ising model coupled to mini- 
mal dynamical triangulations. In analogy to the situation in conventional dynamical 
triangulation we expect that the addition of Ising spins to the graphs leads to the 
possibility of a spin-ordering phase transition at some critical value of the temper- 
ature f3c [10|. Performing simulations on triangulations ranging from 250 to 8000 
vertices we indeed observe all the usual signs of such a transition; a peak in the spin 
susceptibility which narrows and grows with increasing lattice volume, a non-zero 
magnetization for large coupling and a cusp in the specific heat. 

We have attempted to determine the infinite volume critical coupling f3c using 
an extrapolation of the pseudo-critical coupling (3c{A) associated with the peaks 
in the specific heat and the derivative of Binder's cumulant gr- In both cases the 
location of the peaks is expected to approach /5c as |/3c — PdA) \ ~ A~^/'^^^ . This fit 
is made possible by an independent determination of I'dn ~ the height of the peak 
in dgr/dp scales as A^/^'^^ . In Fig. 3 we show the scaling of the pseudo-critical 
couplings together with the infinite volume extrapolation. Our best estimate of the 
critical coupling from this is (3c = 0.2663(3). 

It is interesting to compare this value to the critical temperature of the Ising 
model on other classes of triangulations: (3c ~ 0.2747 (fixed triangular lattice), 
(3c ~ 0.2163 (combinatorial), and (3c ~ 0.1603 (degenerate). It is clear that magnetic 
ordering occurs more readily as the triangulation space is enlarged. 

To obtain the critical exponents of the model we use standard finite size scaling. 
We have already discussed vdn- In the same way a is obtained from the scaling 
of the peak in the specific heat: Cy = cq + ciA°^/'^'^" . From the scaling of the 
magnetization M and the spin susceptibility x the critical point we get (3 and 
7 respectively. The results are shown in Table There we also demonstrate the 
finite size effects present by systematically excluding data for small volume. The 
results are compared to the corresponding values for the Onsager solution of the 
Ising model on a fixed 2D lattice, and the KPZ values for the Ising model coupled 
to 2L'-gravity. It is clear that for the Ising model coupled to the minimal class of 
triangulations the values agree with those of KPZ. An analytic formulation of a 
model with vertices of restricted coordination number has recently been given |12|. 
This model treats squares with coordination number three, four or five. Although 
discrete loop equations can be derived they have yet been solved. 

In this letter we have provided strong evidence that a new class of triangulated 
models, the MDT models, exhibit the same critical behavior as the conventional 
dynamically triangulated random surface models of two-dimensional quantum grav- 
ity. We have demonstrated this both in pure gravity, where the critical properties 
are encoded in the string susceptibility 7^, and in the critical Ising model coupled 
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Figure 3: The infinite volume extrapolation of the pseudo-critical couplings obtained 
from the peaks in the specific heat and the derivative of Binder's cumulant. The data 

points correspond to A = 500, 1000, 2000, 4000 and 8000. For the Binder's cumulant 
fit only the three largest volumes are included in the extrapolation since finite-size 
effects are more pronounced. The rescaling uses the measured value i^dn = 2.9. 



Table 2: The measured critical exponents, udn and a are obtained from the scaling 
of the peaks in dgr/d(3 and Cy respectively. /3 and 7 are measured from scaling at 
Pc, in which case the errors are dominated by the uncertainty in the location of Pc- 
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to gravity where the Ising exponents yield a sensitive test of the critical behavior 
of the model. Our results suggest that putting rather severe cut-offs on the vertex 
coordination number does not affect the universal behavior of the model. This is 
a highly non-trivial observation and implies that it may indeed be possible to find 
physical 2D systems which exhibit this universal behavior []Tl| ]. 

For quantum gravity it implies that large fluctuations in the curvature at the 
scale of the lattice cut-off are unimportant for determining the continuum structure 
of 2D quantum gravity - small curvature defects are sufficient. The pure gravity 
model may be interpreted, similarly to as a gas of strength ±1 curvature defects 
(or vortices) in flat space. In the same way the Ising-type model consists of an 
interacting gas of Ising spins and curvature defects - the gravitational dressing of 
the Ising spins being a function of their interaction with these discrete defects. In 
this context, it may be possible to get some insight into the physical origins of this 
dressing within this simple model. A study of the geometry of spin clusters and its 
relation to the curvature defect distribution is certainly warranted. 
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stimulating discussions. This work was supported in part by the Department of 
Energy U.S. A under contract No. DE-FG02-85ER40237 and from research funds 
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